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Abstract 

We consider the tentative indication of flavor-violating Higgs boson decay h ^ gr recently reported 
in the CMS experiment within the framework of minimal flavor violation. Specifically, we adopt the 
standard model extended with the seesaw mechanism involving right-handed neutrinos plus effective 
dimension-six operators satisfying the minimal flavor violation principle in the lepton sector. We find 
that it is possible to accommodate the CMS h ^ gr signal interpretation provided that the right- 
handed neutrinos couple to the Higgs boson in some nontrivial way. We take into account empirical 
constraints from other lepton-flavor-violating processes and discuss how future searches for the /r —)• ey 
decay and g ^ e conversion in nuclei may further probe the lepton-flavor-violating Higgs couplings. 
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I. INTRODUCTION 


The Higgs boson discovered at the LHC three years ago [|T[ can offer a potential window into 
physics beyond the standard model (SM). The existence of new interactions can bring about 
modifications to the standard decay modes of the particle and/or cause it to undergo exotic 
decays |^. As LHC data continues to accumulate with increasing precision, they may reveal 
clues of new physics in the Higgs couplings. 

The latest LHC measurements of the Higgs, h, have started to expose its Yukawa interactions 
with leptons. Particularly, the ATLAS and CMS Collaborations have observed the decay mode 
h —)■ T^T~ and measured its signal strength to be cr/cTg^ = and 0.91 ± 0.28, respec¬ 
tively In contrast, their direct searches for the decay channel h —?■ have so far come 

up with only upper limits on its branching fraction, B{h —)■ < 1.5 x 10“^ and 1.6 x 10“^, 

respectively at 95% confidence level (CL). Overall, these results are still consistent with 

SM expectations. 


There have also been searches for flavor-violating dilepton Higgs decays, which the SM does 
not accommodate. In this regard, CMS recently reported 0 the interesting detection of a slight 
excess of h ^ events with a significance of 2.4cr. If interpreted as a signal, the excess 

implies a branching fraction of B{h —)■ ^t) = B{h —)■ + B{h —)■ = (0.841q;37)%, 

but as a statistical fluctuation it translates into the bound B{h —)■ /ir) < 1.51% at 95% CL [0. 
In view of its low statistical significance, it is too soon to draw a definite conclusion from this 
finding, but it would constitute evidence of new physics if confirmed by future experiments. 


This tantalizing, albeit tentative, hint of lepton flavor violation (LFV) outside the neutrino 
sector has attracted a growing amount of attention, as the detection of such a process would 
serve as a test for many models [p|-pT[| and could have major implications for upcoming Higgs 
measurements Biig. Subsequent to the h —)■ pr announcement by CMS, its signal hypoth¬ 
esis was theoretically examined in the contexts of various scenarios involving enlarged scalar 
sectors |[T3HT6| or nonrenormalizable effective interactions [p!^-[T7[|. 


In this paper, we follow the latter line of approach which relies on effective operators to 
address LFV in Higgs decay. To handle the LFV pattern systematically without getting into 
model details, we adopt the framework of so-called minimal flavor violation (MFV). Motivated 
by the fact that the SM has succeeded in describing the existing data on flavor-changing neutral 
currents and CP violation in the quark sector, the MFV principle presupposes that Yukawa 
couplings are the only sources for the breaking of flavor and CP symmetries [^, However, 


unlike its straightforward implementation for quarks, there is no unique way to extend the notion 
of MFV to leptons, as the minimal version of the SM by itself, without right-handed neutrinos or 
extra scalar particles, does not accommodate LFV. In light of the fact that flavor mixing among 
neutrinos has been empirically established it is attractive to formulate leptonic MFV by 
incorporating new ingredients that can explain this observation [^. Thus, here we consider 


the SM expanded with the addition of three heavy right-handed neutrinos as well as effective 
dimension-six operators conforming to the MFV criterion.^ The heavy neutrinos are essential for 
the seesaw mechanism to endow light neutrinos with Majorana masses. 


Various scenarios of leptonic MFV have been discussed in the literature |21|-P5| 
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In the next section, after briefly reviewing the MFV framework, we introdnce the effective 
dimension-six operators that can give rise to LFV in Higgs decay, only one of which is relevant 
to h —/iT. In Section we explore the parameter space associated with this operator which 
can yield B{h —)■ /rr) ~ 1%, as CMS may have discovered. At the same time, we take into 
account various experimental restrictions on the Higgs couplings proceeding from the operator. 
Specifically, we impose constraints inferred from the LHC measurements described above as well 
as from the existing data on transitions with LFV that have long been the subject of intensive 
quests, such as /r —)■ ey. We present several sample points from the viable parameter space 
that can account for the CMS’ h —)■ /rr signal interpretation. We also discuss how future 
searches for /i —)■ ey and nuclear /r —)■ e conversion may offer further tests on the interactions 
of interest. Finally, we look at a few other processes that can be induced by the same operator. 
Especially, we hnd that the Z-boson decay Z —)■ /rr can have a branching ratio that is below its 
current empirical limit by merely less than an order of magnitude. We make our conclusions in 
Section m An appendix contains some additional information and formulas. 


II. OPERATORS WITH MINIMAL LEPTON-FLAVOR VIOLATON 


In the SM plus three right-handed Majorana neutrinos, the renormalizable Lagrangian for 
lepton masses can be written as 




0^e)kl ^k,L^l,R^ \i^i)kl^k,R^l,R + H.C. , 


( 1 ) 


where /c, / = 1, 2, 3 are implicitly summed over, W e denote Yukawa coupling matrices, ^ stands 
for left-handed lepton doublets, ^ and represent right-handed neutrinos and charged lep¬ 
tons, respectively, H = iT 2 H* with T 2 being the second Pauli matrix and H the Higgs doublet. 
My is the Majorana mass matrix of Ui and = {i'k,RY, the superscript referring to charge 
conjugation. For the nonzero elements of My taken to be much greater than those of vYy/\/2, 
the seesaw mechanism of type I is operational and generates the light neutrinos’ mass matrix 
rriy = —{v‘^/2)Y^M~^YJ = , where v ~ 246 GeV is the Higgs’s vacuum expec¬ 

tation value, HpMNs denotes the Pontecorvo-Maki-Nakagawa-Sakata (PMNS |^) matrix, and 
rhy = diag(m]^, m 2 , m 3 ) contains the light neutrinos’ eigenmasses. This suggests 


n = 


i \/2 




( 2 ) 


where O in general is a complex 3x3 matrix satisfying 00^ = I, the right-hand side being 
a unit matrix, and can be parameterized as 


0 


O = 




RW = 


S') 

1 

S') 


S') 


S') 


S')' 

2 

S') 


(3) 


with 2 3 r j 2 3 being independent real constants. Hence nonvanishing r )'2 3 dictate how the 
Higgs couples to the right-handed neutrinos in a nontrivial way according to Eq. ( 0 )- Hereafter, we 
concentrate on the possibility that the right-handed neutrinos are degenerate, so that My = All. 
In this particular scenario, only the part of O matters physically 
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The MFV hypothesis [^, ^ then implies that £m is formally invariant under the global 
flavor group Gi = SU(3)l x 0{3)i, x SU(3 )e. This entails that Lk,L, and Ek^n belong to 
the fundamental representations of their respective flavor groups, 

y ^u^R y ^R ^E^R y (4) 

where Vl^e ^ SU(3)i:,^£; and O^, G 0(3)^ is an orthogonal real matrix Furthermore, 

under Gi the Yukawa couplings transform in the spurion sense according to 

i; ^ uuoj, n ^ unvj. (5) 


Due to the symmetry under Gi, we can work in the basis where = \/2 diag(mg, m^, /v 
and the helds i^k,R, and refer to the mass eigenstates. Explicitly, {Ei, E 2 , E^) = (e, /i, r). 
We can then express in relation to Fp 


as 


Lk,L — 


In the standard parametrization 


(h^PMNs)fci D,L 

^k,L 


( 6 ) 




PMNS 


■'512*^23 


-12 "-IS 

“ '-12'523'®13 


'^12 ^23 ^ 12 ^ 23 ^ 13 ^ 


iS 


*-12 *-23 


’12 "-13 
'5i2'523'®13 ® 


^12^23 ^12^23^13® 


iS 

iS 


o o;^3 
'^ 23'-13 
^23 ^13 


diag(e*“i/^e*“^/^l), (7) 


where S and 01^2 are the CP-violating Dirac and Majorana phases, respectively, Cf^i = cos 0f,i, 
and S/^i = sin 6*^;. 

To put together effective Lagrangians beyond the SM with MFV built-in, one inserts products 
of the Yukawa matrices among the pertinent helds to assemble G^-invariant operators that are 

[1. Of interest here are the combinations 


singlet under the SM gauge group |]I9 


A = Y U 

^ ^1/^1/ ^2 ^PMNS *-^PMNS ’ 


B = = diag(l/,^l/J,l/2) , (8) 


where = \/2mj/u. With these matrices, one can generally devise an object A as an inhnite 
power series in them and their products, but it turns out to be resummable into only 17 terms . 
To maximize the new-physics effects, we assume that the right-handed neutrinos’ mass A4. is large 
enough to render the biggest eigenvalue of A equal to unity, which conforms to the perturbativity 
requirement |^. Given that the eigenvalues of B are at most ~ 1 x 10“^, we may 
consequently drop from A all the terms with B, which would otherwise be needed in a study 
concerning GP violation Accordingly, the relevant building block is 


^ — ^i1+^2A + ^4A^ 


(9) 


where in our model-independent approach ^124 free parameters expected to be at most 
of 0{1), one or more of which could be suppressed or vanish, depending on the underlying 
theory. As Im ^^2 4 tiny we can further approximate A^ = A. 
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One could then construct the desired Grinvariant effective Lagrangians that are SM gauge 
singlet. The one pertaining to h —)■ at tree level is given by 


C 


O 


MFV 


(e3) 

RL 


+ H.c,, Og> = 


A2 


( 10 ) 


where the mass scale A characterizes the underlying heavy new-physics and the covariant deriva¬ 
tives VPR = d^H + i[gT^W^ + g'BP)H/2 and V^L = d^L + i{gT^WP — g'Bf’)L/2 contain the 
usual SU(2)i X U(l)y gauge fields and B^ with coupling constants g and g', respectively, and 
Pauli matrices with summation over a = 1, 2, 3 being implicit. There are other dimension-six 


MFV operators involving H and leptons that have been written down 




with being of the form of A in Eq. 




( 11 ) 


and having their own coefficients but these 
operators do not induce tree-level dilepton Higgs couplings. The same thing can be said of the 
comparatively more suppressed i[H^VpH — {VpHy‘^ rr^sEr. In the literature the 
operator HErY^AH^Lr is also often considered {e.g., 0), but it can be shown to be related 
to and the other operators above. Explicitly, employing the equations of motions for SM 
helds [BDl 


one can 


derive 


0g' + H,c. = ^[H^V,H-(V,HyH]{L^Y{^,Y,Y^}L^+iEgYYiAY,E^) 
+ g [HtV,H + (.V,H)^H] [A, YX'] E, 

+ + (E>XYH]ELY[iY,YX]TX 

+ i - 2)mlEX/YH'E, + BXAL,^ 

+ B„n'Aff^//'(j'B'"+jT„»'r)U + H.C.] 


( 12 ) 


plus terms involving quark helds and total derivatives.^ The third and fourth lines of this 
equation, which have [A, , also supply contributions to h —)■ ii', but they correspond to 

small, O {m^ , effects that will be ignored later in Eq. (pAf). 

III. DECAY AMPLITUDES AND NUMERICAL ANALYSIS 


One can express the effective Lagrangian describing the Higgs decays h ^ i for 

£ 7 ^ F as 


Emr = -yu~^PR^' -yi'i^'PR^ + H.c., 


(13) 


The formula for H.c. in the footnote 1 of Ref. has several terms missing and the wrong sign in the 


dipole {(Jpu>) part. These errors have been corrected here in Eq. (|12[) . 


5 















where 3^^^; denote the Yukawa couplings, which are in general complex. Hence the combined 
rate of h ^ £~ , i'~£~^ is 




= r 




+ r 




Svr 


(\y„,? + \yn?)< 


(14) 


where the lepton masses have been neglected compared to The flavor-conserving decay 
h £~£~^ has a rate of / (Svr). 

The MFV Lagrangian in Eq. ([T0|) contributes to both flavor-conserving and -violating Higgs 




EkEi 


iSM 


— A 

— "fcZ -^EkEk 


m 


El' 


2A^v 


K 




V 

(15) 

’ ^ \yEi\ 

for ££' = ep, er, pr and 


are real in our MFV scenario. 

These couplings enter the amplitudes for a variety of lepton-flavor-violating processes, such 
as p —)■ ey, via one- and two-loop diagrams. Therefore, they are subject to the pertinent 
empirical constraints mu, the most stringent of which we list here, assuming that the impact 
of these loop contributions is not much reduced by other new-physics effects. As we sketch in 


Appendix^, the current bound H(/i ey) < 5.7 x 10 |2^ translates into 


+ S-lSXrXer + Kx,. < 5.1 X 10-' . (16) 

where = 0.29. From B{t e'y)exp < 3.3 x 10“® |2^, one extracts ^ 


y„ + 'vl7lxy+ |XT < 5.2x10-*, (17) 

where = 0.03. In these inequalities, we have put more than two different couplings together, 
as they are generally affected by £mfv at the same time, and dropped smaller terms. The 
aforementioned CMS h —)■ pr result under the no-signal assumption implies |]^ 


dlXT+IXT < 3 6 X W . 


(18) 


which is ^ 4 times stronger than the restraint 
and encompasses the range 
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inferred from B{t —)■ pyjexp < 4-4 x 10 ® 
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2.0 X 10-3 < + \yf,r\^ < 3.3 X 10-3 (19) 

implied by B{h —)■ jjLr) = (0.841q; 37)% in the CMS signal hypothesis 0]. 

The information on h —)■ t~ recently acquired by ATLAS [^, and CMS i, |] is 

also useful for restricting new physics in y^^^rr- From the data described in Section |, we may 
require 

\yjyTA < 6-5. 0.7 < ixxxrr < i.s. ( 20 ) 
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-4 


= 894 eV 


31 


and = 7.19 X 10“^ in the SM from the rates 
for = 125.1 GeV. These numbers allow one to see from Eqs 
represent the 2 -loop effects 


ri|, that the 2 -loop contribution to 


where = 4.24 x 10 

and = 257 ke 

and (p!7D, where and 

p —)■ 67 is dominant in constraining whereas the 1- and 2-loop effects on r —)■ 67 are 

roughly comparable. 

We now attempt to attain \y^^\ ~ 0.003 corresponding to the CMS hint of h —?■ pr by 
scanning the coefficients ^124 in ^ + ^ 2 ^ + ^ 4 ^^ which enter the Yukawa couplings 

according to Eq. (P!5[ ) and consequently are subject to the restrictions in Eqs. (0)-(0) and (0)- 




in 


Given that in our MFV scenario Wr oc m^, if £ 7 ^ from this point on we neglect 
comparison to ye^i,eT,^iT, respectively. 

Since A in Eq. (|[) can be realized in many different ways, we consider Erst the possibility that 
the orthogonal O matrix is real, in which case 


A = YY] = 


2M 


U rh r/f 

'-^PMNS '-^PMNS 


( 21 ) 


and the right-handed neutrinos’ Yukawa coupling matrix in Eq. (El) simplihes to Y^, oc rh 


1/2 


somewhat similar to its Dirac-neutrino counterpart Although Gpmns has dependence on the 
Majorana phases as in Eq. (^, they drop out of Eq. (^T]) . 

To proceed numerically, we employ the central values of neutrino mixing parameters from 
a recent fit to global neutrino data . Most of the numbers depend on whether light neutrino 


masses have a normal hierarchy (NH), < m 2 < m 3 , or an inverted one (IH), m 3 < m^ < m 2 . 

Since experimental information on the absolute scale of m 3 2 3 is still far from precise |^, for 
dehniteness we select m^ = 0 (m 3 = 0) in the NH (IH) case. 

With the preceding choices, after exploring the parameter space, we find that \y^T-\ can 
only reach somewhere in the range of (1-2) x 10“^. This is caused by the constraint in Eq. (PT|), 
without which the upper bound < 0.0036 could be easily saturated. Thus, to reproduce 
the signal range in Eq. (|T^, the form of A in Eq. (^Tf ) is not sufficient, and we instead need one 
with a less simple structure, to which we pay our attention next.^ 

A more promising possibility is that the O matrix in Eq. (|^) is complex, which leads to 


A = Y^YJ = — Af Gp^^gmy^oot^i/s^t 
V 


( 22 ) 


As mentioned in the previous section, one can express O = with real antisymmetric 

matrices R and R' defined in Eq. (^. Accordingly, we have 


00^ = e 




2R' 


sinh^r 


r = Yrl + rl + rl , 


(23) 


and so nonzero ri_ 2,3 can serve as extra free parameters that may allow us to achieve the desired 
size of This can indeed be realized, as illustrated by the examples collected in Table |. 


^ A similar conclusion was drawn in Ref. |p^ from a semi-quantitative investigation focusing on an MFV contri¬ 
bution that corresponds to the G term in our study. 


7 


























ai 

TT 

0^2 

TT 

’’1 

’’2 

^"3 

10® 

(GeV-2) 

(GeV-2) (GeV-2) 

ee 

iSM 

ee 

y 

-^^SM 

y 

w»' -7-7- 

TT 

\ye^ 

I3^erl 

\y^.r\ 

10“® 

10"^ 

10“^ 


0 

0 

0.81 

-1.7 

-0.89 

-6.3 

6.2 

5.4 

1.5 

1.2 

0.89 

1.7 

0.3 

3.1 

NH 

0 

0 

-0.86 

1.8 

-0.92 

-7.1 

8.7 

4.5 

1.6 

1.2 

0.87 

2.0 

0.4 

3.5 


0 

0.23 

0.74 

-0.80 

-0.20 

4.9 

-6.7 

-5.9 

0.63 

0.93 

1.3 

1.7 

2.2 

3.2 


0 

0 

0.04 

0.63 

-0.93 

-7.9 

8.8 

2.6 

1.5 

1.2 

1.1 

2.1 

2.8 

3.2 

IH 

0 

0 

0.02 

-0.75 

1.1 

-5.7 

3.8 

8.1 

1.4 

1.1 

0.90 

2.4 

1.3 

3.3 


0.79 

1.3 

-0.61 

-0.79 

1.4 

-5.3 

5.0 

7.6 

1.4 

1.0 

0.84 

1.2 

0.4 

3.5 


TABLE I: Higgs-lepton Yukawa couplings corresponding to sample values of the Majorana phases 01 ^ 2 , 
the parameters ^1^2, 3 of the complex O matrix, and the coefficients 2 4 MFV building block A 

which can yield \y^^ \ ^3 x 10“^. The calculation of the NH (IH) results also relies on the measured 
neutrino mixing parameters in the case of normal (inverted) hierarchy of neutrino masses. 


The flavor-violating Yukawa couplings quoted in the last three columns have followed from their 
dependence on the elements of A determined using the listed sets of ai 2 i ’" 1 , 2 , 3 , and ■^ 1 , 2 , 4 /’^^ 
numbers, along with the central values of neutrino mixing parameters from Ref. 
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again 


with = 0 (m 3 = 0 ) if the light neutrino masses have a normal (inverted) hierarchy. The 
table includes a couple of instances with nonvanishing Majorana phases 032 , which are not yet 
measured and affect A, as in Eq. ([2^) is not diagonal. 


In the table, we also collect the corresponding flavor-conserving Yukawa couplings divided by 
their SM predictions, including for completeness, with y^^ = mjv = 2.08 x 10“®. It is 
obvious that can be altered sizeably with respect to their SM values. Therefore, measurements 
of h —)■ /i’*'/!”, T^T~ with improved precision in the future can offer complementary tests on the 
new contributions. 


Based on our numerical exploration, there are a few more remarks we would like to make. 
First, we have noticed that the viable parameter ranges in the NH case are broader than their IH 
counterparts. Second, in many trials we observe that \yeT\ for the hypothetical signal 

regions, as Table | also shows. This pattern has implications that may be checked empirically in 
the future. Third, in the absence of either ^2 'C 4 fhe maximal \y^^\ is somewhat lower than 

that when ^124 contributing, but at least some or all of the signal values in Eq. (|I^ can 

be accommodated. However, if only ^ 2 , ^ 4 , or ^ 2,4 ^re nonzero, cannot exceed ~0.0018. 

Now, the six sample sets of parameter values in Table | produce branching fractions of p —)■ eq 
and r —)■ /iq in the ranges of (1.4-5.4)xl0“^^ and (1.6-2.0) xl0“®, respectively, if other new- 
physics effects are negligible. The former numbers are within only a few times below the present 
bound H(/i whereas the latter are at least a factor of 20 less than B{t —?■ 

They can be regarded as predictions testable by ongoing or future experiments looking for these 
decays if the CMS’ indication of h —?■ /ir is substantiated by upcoming Higgs measurements and 
the signal range in Eq. (pT]), or part of it, persists with increased data. Especially, the planned 
MEG H experiment on /r —)■ eq, with sensitivity expected to reach a few times 10“^'^ after 3 
years of data taking [Q will probe the above predictions for it. 























ai 

TT 

0^2 

TT 


^2 

^3 

10® 

(GeV-2) 

io®«2/Y 

(GeV-=) 

10® 54/A2 
(GeV-2) 

3^ 

^ee 

iSM 

^ ee 

y 

-^iSM 

y 

w»' q-j- 

•A iSM 

rr 

I3^e/ 

10“® 

IWrI 

10"^ 

\y,r\ 

10"3 

NH 

0 

0 

-0.53 

0.73 

-0.40 

6.0 

-0.7 

-9.5 

0.53 

0.79 

1.1 

0.6 

0.2 

2.7 

0 

0.4 

0.68 

-0.80 

-0.15 

-5.4 

-2.3 

12 

1.4 

1.2 

0.93 

0.3 

0.5 

2.6 

IH 

0 

0 

0.0 

-0.73 

1.1 

-4.7 

-1.9 

11 

1.4 

1.1 

0.96 

0.5 

0.1 

2.5 

0.8 

1.3 

-0.60 

-0.81 

1.4 

-6.5 

9.4 

1.1 

1.5 

1.2 

1.0 

0.1 

0.5 

2.9 


TABLE II: The same as Table except the // —)• 67 and h —)• fifi, tt constraints are replaced with 
their projected fntnre experimental limits, as described in the text. 


As it turns out, if the forthcoming search for /i —)■ 67 still comes up empty, there could yet 
remain viable, but narrower, signal parameter regions. We illustrate this in Table |T|, assuming 
a possible future limit of B{n —)■ 07 ) < 5 x |Q, which amounts to replacing the right-hand 

side of Eq. (|1^) with 1.5 x 10“^, and also imposing the ratios 0.5 < < 1.5 and 

0.8 < Th^rf/'^^rf < 1-2 based on LHC Run-2 projections |^. Since the examples in Table |I| 
yield H(/i —)■ 07 ) = (1.2-4.4) x 10“^^, they may be out of reach of MEG II, and so to probe them 
one will likely need to rely on experiments looking for nuclear p —)■ 0 conversion, which promise 
a greater degree of sensitivity in the long run . As discussed in Appendix |^, the existing 
data on /i —)■ 0 conversion in nuclei are not yet competitive to the current measured bound on 
/i —)■ 07 in constraining the Yukawa couplings. However, we also point out in the appendix that 
planned searches for /i —)■ 0 conversion, such as Mu2E and COMET |^, can be expected to test 
very well the parameter space represented by the examples in Tables | and E 

Finally, we discuss the contributions of £mfv in Eq- m to some other processes. Expanding 
the operator, we have 


^(e3) _ 


B^Pl ( 9nE, - leA^E, + tg^Z^E, + ^ W-v, ] !fh 


\i9mE, 


E,Pl 


iZW„E, iW-d^Ui eA-Z grZ^ 9 

- E-L -^-i +--L T IVAW~ E, 

2 c„ V2 2 c„ 2 c^ + 2 '^'^ ' ' 


(h + v), 

(24) 


where = 9{sl - and c„ = a/ 1 “ -^w = 9v/{2mz) = my^/mz- Evidently, £mfv 

not only induces the already addressed h —>■ W couplings, but also contributes to the two- 
body decays of the weak bosons, Z —)■ ll' and W —)■ ri/;, as well as to three- and four-body 
modes, such as h ^ ££' 7 , , lI'^Z. Since the latter are more suppressed by phase space, we 

deal with only the two-body Z and W decays. The other operators in Eq. ([TT|) can also affect 
Z —>■ ll' and W —)■ rz/^, but here we entertain the possibility that their impact is comparatively 
unimportant. Accordingly, from Eq. (|2^) we derive 


■^Z^EkEi — "^Et, 


/\ / \ 
^ki iz {9lPl + 9rPr) S ’^2y^ ^z ' Pei ~ "^EiPr ^z ' Pe^j 


^Ei > 


M 


W^TU, 


= 


^31 9 


V2A. 




%/2 


AT 


(25) 
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where g^ = gsl^/c^ and we have included the SM terms in these amplitudes. Hence, neglecting 
lepton masses compared to we arrive at 


Z—VjlG 


= r 


Ainm, 


rrir 


Z — VfJjG. 


192 


dSvrm^ 


(26) 


and similarly for Z eT,HT. Thus, for, say, [D^e^l ~ x 10-«. = 2.8 X 10 and 

13^^.,-1 = 0.0032 from Table |, we get 

B{Z eV^) = 6.0 X 10-^^ B{Z e^r^) = 1.1 x 10"®, B{Z = 1.4 x 10"®. (27) 


For comparison, the experimental limits are 


B{Z < 1.7 X 10"® , B{Z < 9.8 x 10“® , 

B{Z < 1.2 X 10"® 


(28) 


at 95% CL. We see that the predicted B{Z —)■ pr) is below its experimental bound by only less 
than a factor of 10. Therefore, Z —)■ pr is potentially more testable than Z —)■ e/i, er, and the 
quest for it can provide a complementary check on £mfv. 

Neglecting lepton masses compared to we also obtain from Eq. ( pS]) 


Z — 


W ^TU 



m 


w 


48e 


2 , 


2 I 

^ 2A4n2 



+ |A 32 p)m^m' 


w 


96A^ev^ 


(29) 


where in the W —>■ te formula we have summed over the 3 neutrino flavors. For the parameter 
values in Table ||, the nonstandard terms in F z^E^Bk being smaller than the 

SM parts by more than 4 orders of magnitude. 


Before ending this section, we would like to note that all the preceding analysis can be repeated 
within the context of the type-III seesaw model with MFV, which is very similar to the type-I 
case addressed in this study if the triplet leptons in the former are as heavy as the right-handed 

E3- However, in the type-H seesaw model with MFV, the Yukawa 


neutrinos in the latter 

coupling matrix of the triplet scalars does not possess the special feature that Y^, has with regard 
to the O matrix that allows to become large enough to explain the CMS h ^ fir signal 
hypothesis. 


IV. CONCLUSIONS 

We have explored the possibility that the slight excess of h ^ jir events recently detected 
in the CMS experiment has a new-physics origin. Adopting in particular the effective theory 
framework of MFV, we consider the SM extended with the type-I seesaw mechanism and an 
effective dimension-six operator responsible for the flavor-violating dilepton Higgs decay. We 
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demonstrate that to account for the tentative h ^ fir signal, with a branching fraction of 
order 1%, the Yukawa coupling matrix of the right-handed neutrinos needs to have a nontrivial 
structure because of the stringent empirical constraints. To illustrate this, we present several 
benchmark points that have survived the restrictions from the existing /i —)■ 07 , r —)■ 07 , 
and h —)■ /r/i, rr data. The viable parameter space can be probed further by upcoming LHC 
measurements and future quests for charged-lepton-flavor violation. Lastly, we examine a few 
other transitions that arise from the same dimension-six operator, among which Z —)■ /rr can 
have a predicted branching ratio merely less than 10 times below its current empirical limit and 
hence potentially also testable in near-future searches. 
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Appendix A: Constraints from p —>• ey decay and p —>• e conversion 


The effective Lagrangian for /i ^ ey can be expressed as 


C 


/i—)-e7 


471-2 


^‘^'“(ClPl + CrPrVF^ 


(Al) 


where a = 1/137 is the hne structure constant, Pl,r = ( 1 T 75 )/ 2 , and Fp^ is the electromagnetic 
held strength tensor. This leads to the decay rate 

A 




am. 


647r^ 


(icip+icHr), 


(A2) 


The Wilson coefhcients Cl,r receive contributions from Higgs-mediated one-loop and two- 


loop O diagrams, 


2 loop 


Given that is real and |A’gg| |T 


hnds [PT| 




one 


q 1 loop 
R 


q 2 loop 
R 


y y 

^ /i/i ^ e/i 

2ml 

0.055 mV 


mh 2 


log^-7 + 


m, 


2m ^ml 


log — - 7 

m^ 4 


efi 


m^ml 


(A3) 


and obtainable from (jj^i°°P' 2 ioop replacements y^^, —)■ yp^. Here we suppose 

that there are no other new-physics contributions that can bring about destructive interfer¬ 
ence with these coefhcients. Thus, putting together these formulas with the latest experimental 
bound @ H(/i —)■ e 7 )g < 5.7 x 10“^^, we arrive at Eq. ([IBD for = 125.1 GeV, which is 


consistent with the most recent measurement 39 


The ehective Lagrangian for /i —)> e conversion in nuclei is 

fanm,, , . 1 „ x 

ea^ ^{Qls^r + > 


r = 


47r^ 


(A4) 
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where q runs over all quark flavors, we have displayed only the most important terms for our 
purposes, and, if yw are the only LFV sources, are already written down in the preceding 


paragraph and [pfl 


9ls 


—2m y* 

O' ^ /ie 

m\v 


9rs 


-2i7i y 

^ efi 

m\v 


The p —)■ e conversion rate in nucleus J\f is then given by 


m. 


yN) _ sr 9 ls f ( jv )„ 

~ ^ mj" 


f A 

‘^capt 

—2m 


an ClDj^ ^(p) ^(p) -in) qin) 
9ls 9lS 


Stt^ 


(A5) 


+ (T —y R) , (A6) 


mN = 


Ny^e 


m\v 


Y f(^) fW = 

^ Jo ? J a 


{N\m„qq\N) 


m 


N = p,n, (A7) 


N 


where Dj^ and are dimensionless integrals representing the overlap of electron and muon 

wave functions for J\f and oj^pt is the rate of muon capture in J\f. Based on the current ex¬ 


perimental limits on p —)■ e transition in various nuclei M and the corresponding overlap 


values |OT, one expects that the J\f = An and Ti data may supply the most 


integral and 

consequential restrictions. The evaluation of —)■ eJ\f) for these two nuclei, respectively, 

requires Dau = 0.189, Dti = 0.087, = 0.0614, = 0.0918, 4^ = 0.0368, 4"^ = 0.0435, 

= 13.07 X 10®/s, and = 2.59 x 10®/s [^, as well as the latest determination of the 
sum of the nucleon matrix elements, = 0.305 ± 0.009 which lies around the lower 

end of the ranges from some of earlier estimates [Q . 

If we impose the measured bound ;B(/iAu —)■ eAu)exp < 7 x 10 “^^ || 2 ^, instead of Eq. (|T6|) , 
but still apply Eqs. ( 0 ), and (13) ! up with l^g^l < 1.6 x 10 which is compatible 

with the finding of Ref. [^. If we use A/” = Ti with B{p,Ti —)■ eTi)exp < 6.1 x 10“^^ [^ , 
instead of Af = An, we get the somewhat stricter \y^p^\ < 1.3 x 10“®. These limitations are 
roughly 5 to 13 times higher than the range of results |3^gp| = (1.2-2.4) x 10“® quoted in Table |, 
demonstrating that the present data on nuclear /i —)■ e conversion are not yet competitive to 
B{p —)> e 7 )exp in restricting especially y , which is also known in the literature |TT], ^ 


Nevertheless, the leading planned searches for p ^ e conversion, Mu2E and COMET, which 
utilize aluminum as the target material [Q, will likely be able to probe the parameter space 


represented by the examples in both Tables | and More precisely, from the sets of sample 
numbers in these tables, together with the aluminum parameters Dai = 0.0362, = 0.0155, 

4f = 0.0167, and 4^ = 0.7054 x lOVs 


we obtain B{fiAl —)■ eAl) = (0.1-9.0) x 10 


-15 


which are within reach of Mu2E and COMET, expected to have sensitivity levels under 10 or 


better after several years of running [34 
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